#### EXERCCIO 9
#### OLIDAN POCIUS
#### Exerccio Programar 1.: QUE FRIO!
celsius<-function(x,round=TRUE)  		
	{
	if (round==TRUE)
	{
	resultado <- round (5/9*(x-32),1)
	}
	  else
		{
		resultado <-  5/9*(x-32)
		}
		return(resultado)
  }

celsius(c(50,100,200))
celsius(c(50,100,200),round=FALSE)

#### Exerccio Programar 2. : Anlise exploratria simultnea de duas variveis
EXPL <- function (x, y, nomes) {
dados<-data.frame(x,y)
names(dados)=nomes
suma<-summary(dados)
correl<-cor(dados)
retorna<-list(suma,correl)
names(retorna)<-c("Sumrio", "Correlacao de Pearson")
plot (y~x,xlab=names(dados[1]),ylab=names(dados[2]),main=paste(names(dados[2]),"versus", names(dados[1])))
linha<-lm(y~x)
abline(linha,col="red")
return(retorna)
}
dados1=c(12,11,34,24,55,56,78,55,15,67,21)
dados2=c(dados1^1.6)
nomes=c("Abacaxi","Melancia")
EXPL (dados1,dados2,nomes)
# $`Sumrio`
#     Abacaxi         Melancia      
#  Min.   :11.00   Min.   :  46.37  
#  1st Qu.:18.00   1st Qu.: 103.32  
#  Median :34.00   Median : 282.08  
#  Mean   :38.91   Mean   : 408.61  
#  3rd Qu.:55.50   3rd Qu.: 617.85  
#  Max.   :78.00   Max.   :1065.00  
# 
# $`Correlacao de Pearson`
#            Abacaxi  Melancia
# Abacaxi  1.0000000 0.9925198
# Melancia 0.9925198 1.0000000

#### Programar 3: ndices de Diversidade de Espcies

diversity <- function(x, modelo="shannon") {
 ab=apply(x,2,sum,na.rm=TRUE)
 pi=t(x)/ab
 pi[pi==0]<-NA
 pi2=pi^2
 logpi=pi*log(pi)
 	if(modelo=="simpson")
 	{
 	diversity=apply(pi2,1,sum,na.rm=TRUE)
 	}
 	if(modelo=="shannon")
 	{
 	diversity=(-1)*apply(logpi,1,sum,na.rm=TRUE)
 	}
 return(diversity)
 }
 

